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Abstract 


Water transport phenomenon in PEM and the mechanism of occurrence and development of a two-phase countercurrent flow with corresponding 
transport phenomenon in the PEM are analyzed. A one-dimensional steady state model of heat and mass transfer in porous media system with 
internal volumetric ohmic heating is developed and simulated numerically. The results show that two dimensionless parameters D and N, which 
reflect the liquid water flow rate and inner heat source in the PEM, respectively, are the most important factors for the water fraction and thermal 
balance in the PEM. The saturation profiles within the two-phase region at various operating modes are obtained. Smaller mass flow rate of liquid 
water and high current density are the major contributions to the membrane dehydration. 
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1. Introduction 


The direct methanol fuel cell (DMFC), which allows direct 
use of an aqueous methanol solution as the fuel and air (O2) 
as the oxidant, is a new kind of proton exchange membrane 
fuel cell (PEMFC). Due to its much simpler peripheral units 
than those of traditional hydrogen PEMFC, which is difficult 
to safely feed the hydrogen in reality utilization, it has been 
attached more and more attention [1]. Based on the experimen- 
tal tests of a DMFC, water and thermal management is one of 
the key factors, which ensure its highly efficient steady opera- 
tion and extend its life-span. There has already been quite a few 
investigations on water and thermal management for a DMFC 
[2-11], and the prevention of dehydration phenomena in PEM 
has also been investigated. However, to our knowledge, there 
is few published literature on the mechanisms of dry-out and 
burnout phenomena in PEM, which take place in the case of a 
certain high current density. Usually, water in a working DMFC 
comes from the aqueous methanol solution fed into the anode, 
humidified air (O2) fed into the cathode and the product of cath- 
ode electrochemical reaction. With the more abundant sources of 
water supply, in the most normal operation conditions a DMFC 
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can work under the higher current density without the threat 
of dehydration of PEM, but in some formidable conditions the 
burnout of PEM was observed in the tests. Hence, to prevent 
dry-out and burnout of PEM under high current density will do 
great help for enhancing the performance and longevity of a 
DMEC. When current density increases to a certain high value, 
evaporation and condensation processes of water may go syn- 
chronously with its transport in a DMFC. Consequently, to find 
out the mechanisms of heat and mass transfer in a DMFC is 
the hard core of water and thermal management. Springer et al. 
[12] and Springer and Wilson [13] discussed the phenomenon 
of two-phase transport in traditional hydrogen PEM fuel cells. 
However, there is still a lack of understanding on the mecha- 
nism of occurrence and development of a two-phase flow in a 
DMEC. In recent years, Wang [14] and his group applied the 
theory of multiphase, multi-component transport in capillary 
porous media into the modeling of the PEMFC, and numerically 
simulated the two-phase flow in porous electrodes considering 
the evaporation and condensation of water. Then some abnor- 
mal phenomena with the concomitant phase change of water, 
such as cathode flooding, dehydration and dry-out of PEM, etc., 
which appear in some formidable operation modes, started to 
be investigated. Since the PEM is very thin and there is only 
a small range of operating conditions between drying out and 
flooding of the cathode [11], it is so hard to carry out the visu- 
alized experiment of water transportation inside PEM. Thus, to 
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Nomenclature 


C concentration of a species 
F Faraday’s constant 

hfg latent heat of evaporation 
I current density 

m mass flux 

p pressure 

q Joule heat flux 

S saturation 

T temperature 

ud permeability velocity 

ule electro-osmotic velocity 
x abscissa 

Z charge number of a species 


Greek symbols 

kinematic viscosity ratio 
non-dimensional length 
porosity 

non-dimensional temperature; contact angle 
permeability 

dynamic viscosity 

kinematic viscosity 

density 

vapor-—liquid interfacial tension 
electric potential 


Sad TEA DV OD 


Subscripts 

ac interface between anode and catalyst layer 
l liquid 

s saturation 

v vapor 


understand the transport process of water in PEM, it is essential 
to develop a theoretical model to analyze the characteristics of 
water transport in PEM, and explore the mechanism of dry-out 
of PEM. In this paper, a one-dimensional steady state model 
of heat and mass transfer in porous media system with internal 
volumetric ohmic heating is developed and simulated and the 
mechanisms of occurrence and development of two-phase zone 
in PEM of a DMFC as a porous media are explored. 


2. The flow and phase change in PEM 


The standard membrane material used in PEM is a fully flu- 
orinated Telfon-base material. Water transport in hydrophilic 
membrane is molecular diffusion. Water molecules exist in the 
pore space of membrane structure. Thus, the water molecular 
transport in hydrophilic membrane can be considered as a diffu- 
sion and permeable process in a porous media with nanometer 
scale. Therefore, the PEM can be regarded as a porous media 
system. To analyze the mechanism of water transport in PEM, 
the PEM porous structure can be simplified as a tube bundle 
made of solid matrix and void channels, as shown in Fig. 1. PEM 
is heated from its cathode boundary by heat generated in cath- 


a b 
(a) ; (b) 
Dehydration 
Cc B Z 
PEM Se 
(c) (d) 
Dehydration PEM 
PEM 
Zone Dryout 
PEM Spreading 


Fig. 1. Evolution of two-phase zone and dryout in PEM (a) Nucleation of water 
bubble, (b) The formation of a two-phase zone, (c) Spreading of two-phase zone 
and (d) Dry-out of PEM. 


ode electrochemical reaction, while ohmic heating contributes 
to its internal heat source. Under the normal operating condi- 
tions, void channels of porous PEM are filled with liquid water 
and liquid water in the void channel C will be dragged along the 
channel wall toward cathode side by capillary force. Owing to 
the current density distribution in cross-section of PEM is usu- 
ally non-uniform [15], therefore, when the local current density 
in a certain zone becomes excessively high, the Joule heat pro- 
duced will be high enough to activate the nucleation of water and 
then a bubble zone B will occur as shown in Fig. la. The water 
transport in channel makes this evaporation process of water to 
be continued. Under the function of local Joule heating, the bub- 
bles in the channel will gradually grow up. The whole process 
seems to be that vapor bubbles move in the opposite direction 
with liquid water flow in C, while the liquid water in front of the 
bubble is dragged towards cathode side along the channel wall, 
and thus form a two-phase countercurrent flow in porous PEM 
[16]. At the same time, the reduction of liquid water may cause 
the decrease of local electrical conductivity, which improves the 
local Joule heating. The higher Joule heating can increase the 
superheat of water in the cathode side of void channels C and 
cause a new vapor flow. Thus, the two-phase coexistent zone in 
porous PEM will appear. Meanwhile, vapor flow is continuously 
cooled by contact with liquid water in its front and condenses 
gradually at the interface of liquid—vapor. For a given current 
density, vapor flow will stop at a certain place of PEM due to 
the equilibrium between evaporation and condensation. Thus, a 
local transient two-phase dehydration zone is appeared in PEM, 
as shown in Fig. 1(b). 

With the increase of the local current density, more vapor bub- 
bles generate and the dehydration zone will spread unboundedly 
in PEM. On the other hand, the flux of the liquid water dragged 
by electro-osmotic force increases with the increase of local cur- 
rent density, which may enhance the cooling effect on vapor flow. 
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Fig. 2. A burnout sample of PEM. 


Consequently, the development of local two-phase zone with the 
increase of the current density depends on these two opposite 
effects as Joule heating and electro-osmotic drag, respectively. 
The local dehydration zone may spread with the increasing cur- 
rent density, as shown in Fig. 1(c). Under the most of working 
modes at high current density the normal operation of DMFC 
can still be conducted. The dehydration zone will not cover the 
whole cross-section of the interface between PEM and cathode 
catalyst layer until the local current density reaches to a certain 
high value. At this moment, most of the void channels in PEM 
are filled with vapor and PEM is seriously dehydrated owing to 
liquid water flow collapsing in its cathode side. It causes a large 
contact resistance and leads to a quick increase of local Joule 
heat flux. Finally, the temperature of the solid matrix quickly 
increases, which is the so-called dry-out of PEM, as shown in 
Fig. 1(d). The dry-out of PEM not only disconnects the inter- 
nal circuit of a DMFC and causes burnout of PEM as shown in 
Fig. 2. 


3. Analytical model and simulation 


The cathode catalyst layer, PEM and anode catalyst layer 
combine a porous media system. According to the working 
condition of a DMFC, the system is subjected to a tempera- 
ture difference such that the temperature at the cathode catalyst 
layer is significantly higher than the saturation temperature of 
the interstitial fluid while the anode catalyst layer is maintained 
below the saturation temperature T; as shown in Fig. 3. 

In this porous media system under the boundary and internal 
heating, it has been shown that three distinct regions will appear; 
that is, vapor zone, two-phase zone and liquid zone [17—19]. The 
thickness of PEM is l=X3 — X; and Xe denotes the start point 
of two-phase zone. 


3.1. Water fraction distribution in PEM 


To simplify the analysis some assumptions are made as fol- 
lows: 


Anode catalyst layer PEM Cathode catalyst layer 


Liquid zone 1 Two-phase zone ! Vapor zone 


x, Xp (X,) X% X x 


Fig. 3. Temperature distribution in PEM. 


(1) The local thermodynamically equilibrium in porous PEM is 
maintained. 

(2) Catalyst layers are very thin; they can be treated as the 
boundaries of PEM. 

(3) The flow rate of fluid in porous PEM is low, the inertia force 
could be neglected and the effect of gravity is also neglected. 

(4) The transport properties of liquid and vapor are constant. 
Liquid water is incompressible and vapor is regarded as 
ideal gas. 


The water velocity from anode to cathode at the liquid zone 

in PEM can be calculated by Schlögl equation as: 
ko dọ kpdpı 

am AGIT i akr i + ula (1) 
where the first term of right side is the velocity caused by electro- 
osmotic force and the second is Darcy velocity under a pressure 
gradient. F is Faraday constant, C; is the iron density, Zş is the 
charge number of iron for a species, ¢ is the electric potential, 
and kg, kp are the electro-osmotic coefficient and permeability 
of the water in PEM, respectively. 

The total water velocity in PEM can be expressed as 


uj = Ud + Ue (2) 
The mass flow rate of the liquid water can be calculated as 
mı = pill (3) 


The continuous equation is 


dom) _ 9 


i (4) 


The energy equation under the local thermodynamically equi- 
librium is: 
T d 


d 
p\Cpiui—— 


dT 
oS qa etg +4 (5) 


where Cp is the specific heat capacitance of water, keff is the 
effective thermal conductivity of PEM and q is inner volumetric 
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heat flux produced by Joule heating, 


q=i— (6) 


The relationship between the current density and the electri- 
cal potential as well as electrical conductivity is 


On ge + FC pu T) 


i= 
in which the electrical conductivity om can be expressed as 
Om = F’ Dp T as (8) 


The current continuous equation is 


Te 9 
a (9) 

Liquid water in PEM is mainly existed by the form of Hydrate 
proton (Hj 0), the diffusion of water in PEM can be neglected. 
Thus, the mass flow rate of liquid water can be calculated by 
piu. Assuming xc is the coordinate of the interface between 
liquid zone and two-phase zone. Define the water fraction a; at 
the location x in PEM as: 

x 


(10) 
m 
where m is total mass flow rate of water and mý is liquid flow 
rate at the location x in PEM. 
The average velocity caused by electro-osmotic force and 
average Darcy velocity in PEM are 
Ag 


ie = OZ 6C pF (11) 
üle = — =; 
e ae 


k = 
ing = 2 PR (12) 


Hl l 


where ps is the pressure at x,. Thus, the total average liquid water 
velocity in PEM is 


it) = ita + te (13) 


The following non-dimensional parameters are defined as 


x 
Non-dimensional length L= 7 
; : ; uj — üle 
Non-dimensional velocity U=- — 
uj — Ule 
: f m 
Non-dimensional flow rate M = — 
pıtı 
i : T— Ts 
Non-dimensional temperature 0= 
Tac Ts 
Non-dimensional pressure P= ARE 
Pac — Ps 


1 50 
0.9 45 
0.8 40 
0.7 35 
0.6 30 

= 

P05 25 
0.4 20 
o3] 15 
0.2 \ 10 
0.14 5 
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Fig. 4. The start point length of two-phase zone (D = 1). 


By use of above non-dimensional parameters, Eqs. (4) and 
(5) become 


VU = —=0 14 
JZ (14) 
de opinya (15) 
dL? dL 


where D= Cpml/Kef and N= qÊ Tac — T;)Keff are two non- 
dimensional parameters which reflect the change of water mass 
flow rate and inner heat source in PEM, respectively. 

The boundary conditions are 


L=0, P=1, 6=1 
(16) 
L=L,, P=0, 6=0 


In this paragraph, the water transport phenomena are empha- 
sized and the main objects are the mechanism and development 
of a two-phase countercurrent flow as well as water fraction dis- 
tribution in PEM. The numerical simulation was conducted [21] 
and some calculated results are as follows. 

The start point of two-phase zone length Le (correspond- 
ing to le =X, — X1) and water fraction distribution in PEM are 
shown in Figs. 4—7. It is shown that the start point appears ear- 


1.5 - 7 7 — 1 T 1 T 60 
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1 440 
o 
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0.5 420 
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N [x(-1)] 


Fig. 5. The start point length of two-phase zone (D = 10). 
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Fig. 6. The water fraction in PEM (D= 1). 


lier at the small water flow rate (D = 1) when electrical current 
increases and it causes water fraction in PEM decreasing fast. 
On the other hand, when the electrical current is small (cor- 
responding to small value of N), the water fraction in PEM 
is satisfactory to normal operation condition and the dry-out 
in PEM is never happened. The decrease of the thickness of 
PEM | will decrease the value of N and it is good for the safe 
operation. 


3.2. The length of two phase zone 
The two-phase flow in porous structure of PEM under pres- 


sure gradients can be determined by a modified Darcy’s equation 
as 


ae (17) 


_ __ w dpy 


18 
uv dx i 


The electro-osmotic velocity of liquid water is 


(19) 


Ja=8.47*10e-3 
— N=-15 


--— N=-50 
---- N = -100 


Fig. 7. The water fraction in PEM (D = 10). 


In the above equations, ky and «Kry are the relative perme- 
ability of the liquid water and vapor, respectively and krọ is the 
relative electro-osmotic coefficient of the liquid water. 

The mass flow rate of the liquid and vapor phases in two-phase 
zone are 


mı = pi(uid + ule) (20) 
My = pyly (21) 
The energy and mass conservation equations in two-phase 
zone are 
vi 
q = Mh g = — (22) 
Om 
m +m, =0 (23) 
Solving Eqs. (17)-(23), the following expressions for vapor 


pressure, liquid water pressure and capillary pressure gradients 
are obtained respectively as 


d 
See (24) 
dx KKryh fe 
d d 
Be ar 25) 
dx KKM fe Krl dx 
dpe q Vy vI Kro do 

= Z fC ¢F— 26 
dx Kh fe \ Krv j Krl Krl py dx (26) 


To analyze the length of two-phase zone in PEM, the term 
of water saturation S, which is commonly used in the heat and 
mass transfer in porous media, is introduced in this paper. In 
porous PEM the water saturation S is defined as the volumetric 
percentage of liquid water in a pore volume of PEM with the 
thickness of X — X1. The relationship between water saturation 
S and capillary pressure p, in porous media can be correlated by 
the following expression [22] as 


Pc an (27) 


o coSs he \E 


fS) = 


To simplify the solution, the following non-dimensional 
parameters are introduced: 


. ; xZsCrFAġ K\ 1/2 
Non-dimensional length 6 = —~——— (=) 
o COS be E 
Non-dimensional viscosity and permeability p= 2 
Vy 
2 
Ky] 
= 1/2 
Non-dimensional pressure Py = Pv ~ Ps (=) , 
o cos be \E 
= Pi — Ps (K\ 1/2 
n= Poh (8 
o cos be \é 
: : qvy 
Non-dimensional heat flux oO = ——— 
Kh fpZ fC FAG 


where ø is surface tension of water and @, is contact angle. 
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By use of Eq. (26), Eq. (27) can be rewritten in the following 
dimensionless form 


7 1 8B ds m 
= o( >+ ) =f 5 (28) 


Krl 


dô 
Solving the Eq. (28), we obtain 


ds = —o((1/kry) + (B/kn)) — y (29) 
dô f! 


in which ky, Kry and f are selected based on Ref. [20] as 


f = 1.41701 — S) — 2.12001 — 5)? + 1.263(1 — 5) (30) 
Kr = s? (31) 
Ky = (1 — SP (32) 


Thus, the two-phase zone length can be obtained by integrat- 
ing Eq. (29) from X=X2, S= 1 to X=X4, S=0 as 


0 f! 
oy = 1 dS 
1 —@(O/kry) + (/ka)) — y 

When the saturation gradient is equal to zero within PEM, 
the development of two-phase zone goes steadily and saturation 
will not be changed, thus, there is no vapor zone appearing in 
the PEM. The dimensionless length of two-phase zone 6; cor- 
responding to the distance of |X4 — X2| can be defined as the 
theoretical length of two-phase zone, while the actual length of 
two-phase zone existing in PEM is |X3 — X2|. Obviously, if the 
theoretical dimensionless length of two-phase zone |X4 — X2| is 
less than the actual length of two-phase zone |X3 — X2|, that is, 
the interface between vapor zone and two phase zone is located 
inside of PEM, the liquid water flow collapse in its cathode side, 
the dry-out will exist and the normal operation of cell will be 
stopped. If the thickness of PEM is far less than the theoretical 
length of two-phase zone ô+, the interface between vapor zone 
and two-phase zone is always located outside of the PEM, the 
dryout of PEM would not occur at all. These results can explain 
the existed phenomena reported in literature well, which show 
that when DMFC works at high current condition, the dry-out 
is easy to appear and when a DMFC operates in a low current 
density or the thinner PEM is used, the dry-out and burnout of 
PEM is hard to occur. 

It is found that the parameter |—w/y| is a key factor to analyze 
the characteristics of two-phase zone. A numerical method was 
proposed to calculate the water saturation S and 4; [21]. Some 
calculated results are as follows. 

The saturation distribution curves of liquid water in theoret- 
ical two-phase zone of PEM under various operating modes are 
shown in Fig. 8. It shows that at any position where distance 
6=X — X2 is less than 0.546,, the saturation increases with the 
increase of |—w/y|, while at other position where ô is larger than 
0.546,, the saturation decreases with the increase of |—@/y|. Due 
toga i’, Adi and |—w/y| x i, the increase of current density 
can cause the simultaneous augment of both Joule heat flux and 
electro-osmotic flux of liquid water. The evaporation of liquid 
water in two-phase zone by Joule heating may reduce satura- 
tion, whereas, the electro-osmotic flux supplies liquid water to 
two-phase zone of PEM and causes the increase of saturation. 


Fig. 8. Liquid saturation profiles with various operating modes. 


Within the area where ô is less than 0.546,, the higher value 
of saturation indicates that the local dehydration of PEM is not 
serious. This means that the increase of current density has a 
little impact on the augment of Joule heating. When current 
density increases, the effect of electro-osmotic flux on saturation 
is strongger than that of Joule heat flux in this area. Conversely, 
within the area where ô is larger than 0.546,, the lower value of 
saturation indicates that the local dehydration of PEM is very 
serious, the increase of current density has a large affectation on 
the augment of Joule heating, so the effect of electro-osmotic 
flux of liquid water on saturation is weaker than that of Joule 
heat flux in this area when current density increases. Therefore, 
creating appropriate boundary conditions in anode and cathode 
electrodes to maintain a proper length of two-phase zone in PEM 
has a positive effect on improving performance of a DMFC, 
which is also an effective way to solve some water management 
problems, such as cathode flooding and dehydration of PEM. 
In addition, the evaporation of liquid water in two-phase zone 
is helpful to prevent burnout of PEM in high temperature by 
drawing out the Joule heat. 


3.3. Pressure distribution in two-phase zone of PEM 


The liquid and vapor pressure profiles in two-phase zone of 
PEM with various operating modes are shown in Fig. 9. For 
all cases, it is evident that the largest vapor and liquid pressure 
gradients occur near the liquid zone and vapor zone boundaries, 
respectively. This is due to the low relative permeability near 
the endpoints. It is also noted that the middle parts of liquid and 
vapor pressure profiles approach gradually to be horizontal with 
the increase of |—w/y| (i.e., current density). 

Under the low current density conditions, electro-osmotic 
drag contributes a little to supply liquid water into two-phase 
zone, while the capillary force dominates the process. Owing to 
the small average saturation in low current density, the partial 
pressure gradient of vapor in two-phase zone of PEM is high. In 
the case of high current density, electro-osmotic drag contributes 
a lot to the process and counteracts a part of liquid pressure 
gradient. 
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Fig. 9. Liquid and vapor pressure profiles with various operating modes. 


4. Conclusions 


A one-dimensional porous model for water transport in PEM 
is developed in this work to analyze water transport charac- 
teristics in the PEM. Some major results are obtained and 
summarized below: 


(1) A two-phase zone and vapor-liquid countercurrent flow 
exist in water saturated porous PEM under some operating 
conditions. The flow and heat transfer may prevent burnout 
of PEM under high current density condition. 

(2) The start point of two-phase zone and water fraction distri- 
bution in PEM are analyzed and calculated. The start point 
appears earlier at the small water flow rate and it causes 
water fraction in PEM decreasing fast. When the electrical 
current is small the water fraction in PEM is satisfactory to 
normal operation condition and the dry-out in PEM is never 
happened. 


(3) The theoretical length of two-phase zone of PEM is an 
important non-dimensional parameter to perform the water 
and heat management of DMFC. 

(4) The saturation profiles within the two-phase region at vari- 
ous operating modes are obtained numerically. The results 
can help to analyze and optimize the performance of a 
DMFC operating in its ohmic polarization region. 
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